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Abstract. We investigate the formation of carbon stars 
as a function of the stellar mass and parent metallicity. 
Theoretical modelling is based on an improved scheme 
for treating the third dredge-up in synthetic calculations 
of thermally pulsing asymptotic giant branch (TP-AGB) 
stars. In this approach, the usual criterion (based on a con- 
stant minimum core mass for the occurrence of dredge-up, 
M™ in ) is replaced by one on the minimum temperature at 
the base of the convective envelope, , at the stage 

of the post-flash luminosity maximum. Envelope integra- 
tions then allow determination of M™ ln as a function of 
stellar mass, metallicity, and pulse strength (see Wood 
1981), thus inferring if and when dredge-up first occurs. 
Moreover, the final possible shut down of the process is 
predicted. 

Extensive grids of TP-AGB models were computed 
by Marigo (1998ab) using this scheme. In this paper, we 
present and discuss the calibration of the two dredge- 
up parameters (i.e. efficiency A and T dred ) aimed at re- 
producing the carbon star luminosity function (CSLF) in 
the LMC, using TP-AGB models with original metallicity 
Z = 0.008. In addition to this, the effects of different input 
quantities on the theoretical CSLF are analysed. It turns 
out that the faint tail is almost insensitive to the history of 
star formation rate (SFR) in the parent galaxy, in contrast 
to the bright wing which may be more affected by the de- 
tails of the recent history. Actually, we find that the faint 
end of the CSLF is essentially determined by the temper- 
ature parameter T^ rod . Once the faint end is reproduced, 
the peak location is a stringent calibrator of the efficiency 
parameter A. The best fit to the observed CSLF in the 
LMC is obtained with A = 0.50, logT b drcd = 6.4, and a 
constant SFR up to an age of about 5 x 10 8 yr. This re- 
cent drop of the SFR is invoked to remove a slight excess 
of bright carbon stars otherwise predicted. 

A good fit to the observed CSLF in the SMC is then 
easily derived from the Z = 0.004 models, with a sin- 
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gle choice of parameters A = 0.65, logT^ Icd = 6.4, and a 
constant SFR over the entire significant age interval. The 
result for A is consistent with the theoretical expectation 
that the third dredge-up is more efficient at lower mctal- 
licities. 

Key words: stars: evolution - stars: AGB and post-AGB 
- stars: mass-loss - stars: carbon - Magellanic Clouds 



1. Introduction 

During the TP-AGB phase, recurrent episodes of dredge- 
up can significantly affect the surface chemical composi- 
tion of low- and intermediate-mass stars. These episodes 
are usually known as "the third dredge-up" . The basic 
mechanism of the mixing to the surface can be explained 
on the basis of the standard structure equations applied 
to stellar models in quasi-static evolution (Iben & Ren- 
zini 1983). During the peak of every flash in the helium 
shell, part of the nuclear energy is converted into ther- 
mal energy, and the remaining spent as mechanical work 
for the expansion of the overlaying matter. Material at 
the hydrogen- helium (X — Y) discontinuity and beyond is 
pushed outward so that the H-burning shell is temporarily 
extinguished. At the same time, the base of the convective 
envelope moves inward because of the progressive local 
raising of the radiative temperature gradient over the adi- 
abatic one, triggered by the increase of the energy flux and 
the concomitant decline of the temperature consequent to 
the expansion. 

Dredge-up occurs if the convective envelope can pen- 
etrate into the inter-shell region containing material that 
has just undergone helium burning during the flash. A cer- 
tain amount of inter-shell matter is then brought up to the 
surface, and pollutes the envelope with newly synthesised 
nuclear products. The dredge-up material mainly consists 
of helium, carbon, oxygen and traces of heavy elements 
produced by s-process nucleosynthesis (see Lattanzio & 
Boothroyd 1997 for an extensive review of the topic). 
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The formation of single carbon (C-) stars is usually 
ascribed to the enrichment of 12 C surface abundance pro- 
duced by the third dredge-up. The topic is important for a 
series of reasons. First, C-stars are present in many Local 
Group galaxies, indicating the presence of relatively young 
stellar populations. The precise age and metallicity inter- 
vals relevant for the formation of C-stars, however, are far 
from being fully assessed from both theory and observa- 
tions. Secondly, these stars are crucial for the interpreta- 
tion of a series of chemical abundance data. For instance, 
it is believed that most SiC grains found in primitive me- 
teorites are formed in extended carbon-rich envelopes of 
C-stars (Gallino et al. 1990). The emission features de- 
tected in the infrared spectra from extended circumstcllar 
shells around AGB stars are significantly affected by the 
composition of the dust grains (Ivezic & Elitzur 1995). The 
variety of chemical abundance patterns in planetary neb- 
ulae probably reflects the different efficiency and duration 
of dredge-up episodes in the progenitor stars, as well as 
the possible burning of CNO isotopes at the base of their 
envelopes (the so-called hot-bottom burning or envelope 
burning; see Marigo et al. 1998 and references therein). 

In recent years many investigations have been devoted 
to the analysis of the occurrence and related efficiency of 
the dredge-up, and to the nucleosynthesis driven by ther- 
mal pulses (Boothroyd & Sackmann 1988cd; Frost & Lat- 
tanzio 1996; Straniero et al. 1997; Lattanzio & Boothroyd 
1997; Herwig et al. 1997; and references therein). The 
main challenge of the theory in this case is offered by 
the so-called carbon star mystery, first pointed out by 
Iben (1981). Observations of AGB stars in the Large and 
Small Magellanic Clouds (LMC and SMC) indicate that 
the luminosity range in which carbon stars are found ex- 
tends down to rather faint magnitudes (A/boi <! — 3), 
thus suggesting the evolution from low-mass progenitors 
(Mi ^1 — 3 Mq). In contrast, complete AGB models eas- 
ily explain the formation of carbon stars at much higher 

luminosities (Mboi < —5 6.5), but fail to reproduce 

the faint part of the observed distribution, since dredge- 
up of carbon is hardly found to occur during the TP- AGB 
evolution of low-mass models. 

Although some improvements have recently been at- 
tained (Frost & Lattanzio 1996; Straniero et al. 1997; Her- 
wig et al. 1997) the present understanding of the dredge- 
up process is still unsatisfactory, mostly because of the 
uncertainties in the treatment of the mixing process and 
the definition of convective boundaries. In Sect. || we give 
an overview of the present state of art emerging from the 
available results of complete modelling of the third dredge- 
up, and discuss the role of analytical TP- AGB calculations 
in possibly giving useful theoretical indications. 

In the present study, the treatment of the third dredge- 
up in synthetic TP- AGB calculations has been signifi- 
cantly improved. A simple technique based on integrations 
of static envelope models is used in order to infer if and 
when convective dredge-up is likely to occur during the evo- 



lution of a TP- A GB star. In this way, we abandon the as- 
sumption of a constant parameter M™ m , usually adopted 
in synthetic calculations (Groenewegen & de Jong 1993; 
Marigo et al. 1996a), and the dependence on the current 
stellar mass and metallicity is taken into account. The 
method is described in Sect. ||. 

A primary target of theoretical models of the AGB 
phase is the interpretation and reproduction of the ob- 
served CSLFs in the Magellanic Clouds (Sect. g). This 
is because a sole observable holds a number of profound 
implications for different physical processes occurring in 
AGB stars (e.g. nucleosynthesis, mixing and surface chem- 
ical enrichment, mass- loss), and sets at the same time 
severe constraints both on theoretical prediction of stel- 
lar parameters (e.g. evolutionary rates, luminosities, life- 
times, initial mass of carbon stars' progenitors), and pos- 
sibly on properties of the related stellar populations to 
which carbon stars belong (e.g. the history of star for- 
mation) . Section [| illustrates the derivation of theoretical 
CSLFs. 

In this study the formation of carbon stars is inves- 
tigated with the aid of the parameterized treatment of 
the third dredge-up. In Sect. ^, the dredge-up parameters 
(efficiency A, and base temperature T^ led of the envelope 
at the post-flash luminosity peak) are first determined by 
demanding that the theoretical CSLF, calculated for ini- 
tial metallicity Z — 0.008, suitably fits the distribution 
observed in the LMC. This calibration is important not 
only for deriving theoretical indications on the dredge-up 
process itself, but also because it gives more reliability to 
the analysis on the surface chemical composition of AGB 
stars and related applications (e.g. estimation of chemical 
yields, and chemical abundances in planetary nebulae). 
The sensitivity of the results to the underlying history of 
star formation and the initial mass function (IMF) is also 
explored. 

The same procedure is then applied to find the proper 
pair of parameters to reproduce the CSLF in the SMC 
(Sect. |). In this case evolutionary calculations for initial 
metallicity Z — 0.004 are used. The differences between 
the LMC and SMC distributions are discussed. Section || 
summarises our main conclusions. 

2. The role of synthetic calculations 

The usual analytical treatment of the third dredge-up re- 
quires the knowledge of three basic input quantities: 

— the minimum core mass for dredge- up to occur Af™ m ; 

— the dredge-up efficiency A; 

— the chemical composition of the convective inter-shell. 

A TP- AGB star can experience the third dredge-up when 
its core mass has grown over a critical value Af" lm , as 
indicated by theoretical analyses of thermal pulses (e.g. 
Wood 1981; Iben & Renzini 1983; Boothroyd & Sack- 
mann 1988cd; Lattanzio 1989). Estimates of M c min based 
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Table 1. The third dredge- up in low-mass stars according to full TP-AGB calculations. For given stellar mass Mj, 
envelope abundances (in mass fraction) of hydrogen X and helium Y, and adopted mixing-length parameter a, the 
core mass M™ m at the onset of dredge- up (if found) is given. The value of the bolometric magnitude, Mboi(M — > C), 
corresponds to the pre-flash maximum quiescent luminosity marking the transition to carbon star. The estimates inside 
parenthesis refer to the faintest luminosity during the post-flash luminosity dip. 
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on complete evolutionary calculations are still a matter of 
debate, since this quantity depends on the delicate inter- 
play between model prescriptions (e.g. opacities, nuclear 
reaction rates, convection theory, mass-loss, time and mass 
resolution) and intrinsic properties of the star (e.g. chem- 
ical composition, envelope mass, strength of the thermal 
pulses) . In Table |l| we report some values of M™ n taken 
from different authors. It turns out that the results vary 
remarkably, even among models with the same total mass 
and metallicity. The third dredge-up is often not found in 
low-mass models, and if it is, the corresponding metallic- 
ities are quite low and the transition luminosities [see the 
Mboi(M — > C) entry] are in most cases still brighter than 
the faint wing of the CSLFs in the Magellanic Clouds. 
Nevertheless, a certain trend can be extracted from the 
comparison of the results, i.e. the onset of the third dredge- 
up is favoured at higher stellar mass M, lower metallicity 
Z ', and greater efficiency of envelope convection (obtained 
e.g. by increasing the mixing- length parameter a). 

The parameter A expresses the dredge-up efficiency. It 
is defined as the fraction of the core mass increment over 
the preceding inter-pulse period, AM C , which is dredged- 
up to the surface by the downward penetration of the con- 
vective envelope: 

_ AMd redge 

where AMd re dg C is the mass of the dredged-up material. 



The question of the true value of A is also troublesome. 
The effective penetration of the envelope is crucially af- 
fected by the treatment of the instability of matter against 
convection, still poorly understood. In the standard ap- 
proach, the boundary of a convective region is defined by 
the Schwarzschild criterion, i.e. where the inward grav- 
ity force is balanced by the buoyancy force. However, at 
this point the velocity of the convective eddies is not null, 
so that they can effectively penetrate into formally sta- 
ble regions. This phenomenon is usually known as convec- 
tive overshoot. Over the years, several attempts have been 
made in order to include this effect into evolutionary cal- 
culations (see Chiosi et al. 1992 and references therein). 
However, at present it appears that there is no real agree- 
ment among different authors on the criterion to define 
the maximum inward penetration of the envelope at ther- 
mal pulses during the AGB. Moreover, the results seem to 
depend quite sensitively on technical and numerical de- 
tails (see Frost & Lattanzio 1996 for a discussion of this 
aspect). 

The heterogeneous situation is well represented by the 
two following examples. On the one side, Straniero et al. 
(1997) obtain rather efficient dredge-up just using the 
Schwarzschild criterion and without invoking any extra- 
mixing, but drastically increasing both time and mass res- 
olution during calculations. On the other hand, Herwig 
et al. (1997) find deep extra-mixing during dredge-up on 
the basis of a parameterized description of the overshoot 
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region, built up to mimic the indications from hydrody- 
namical simulations (Freytag et al. 1996) of convection in 
stars of a different kind (i.e. main sequence stars and white 
dwarfs). 

As far as the efficiency of the third dredge-up is con- 
cerned, it is not easy to summarise the results from dif- 
ferent authors. However, suffice it to recall that recur- 
rent estimates, derived from TP-AGB calculations of low- 
mass stars, give at most A ~ 0.3 — 0.4 and often even 
smaller (Iben k Renzini 1984; Lattanzio 1986, 1987, 1989; 
Boothroyd & Sackmann 1988d; Straniero et al. 1997). 

In conclusion, the present state of art is the following. 
On the one hand, most authors still do not find signifi- 
cant dredge-up in their calculations (Vassiliadis & Wood 
1993; Wagenhuber 1996; see also Table 0), or assume the 
efficiency A as an input quantity (Forestini & Charbon- 
nel 1997). On the other hand, the difficulty of a too poor 
dredge-up in low-mass stars seems to be overcome to some 
extent (Boothroyd & Sackmann 1988d; Frost & Lattanzio 
1996; Straniero et al. 1997; Herwig et al. 1997), but exten- 
sive calculations over a wide range of stellar masses and 
metallicities are required to test the results by means of 
a meaningful comparison with observations. Moreover, re- 
sults are quite different even for similar initial conditions 
(e.g. M and Z). 

A powerful investigation tool complementary to full 
calculations is offered by synthetic TP-AGB models, 
which have been notably improved and up-graded in re- 
cent years (Groenewegen & de Jong 1993; Marigo et 
al. 1996a, 1998; Marigo 1998ab). As far as the third 
dredge-up is concerned, synthetic TP-AGB calculations 
clearly suggest that a higher efficiency is required (e.g. 
A ~ 0.65 — 0.75) in order to reproduce observational con- 
straints (e.g. the CSLF in the LMC), a goal still missed 
by complete modelling of the TP-AGB phase. 

The use of the parameters M™ m and A plays a rele- 
vant role in analytical models of the TP-AGB phase. They 
both significantly affect the predicted luminosity functions 
of carbon stars, in particular the position of the peak and 
the low-luminosity tail of the distribution. Lowering M™ m 
favours the formation of carbon stars of lower masses and 
at the same time results in a longer carbon star phase; 
increasing A corresponds to a higher fraction of low mass 
stars which become carbon stars at fainter luminosities. 
In Marigo et al. (1996a) it was simply assumed that both 
M" lm and A were constant, thus neglecting their likely de- 
pendence on the stellar metallicity, total mass, core mass, 
and mixing length parameter. Though being a rough ap- 
proximation, such a choice was meant to provide useful 
theoretical indications. The calibration of the dredge- up 
parameters was based on the goal of reproducing the ob- 
served CSLFs in the LMC. It resulted that the suitable 
values were M c min = 0.58 M and A = 0.65. These val- 
ues are very close to those derived by Groenewegen & de 
Jong (1993). Basing on these results, Forestini & Char- 
bonnel (1997) fixed A = 0.6 in their full calculations, in 



order to estimate the chemical yields from intermediate- 
mass stars. 



3. The minimum core mass for dredge-up 

In this section we illustrate the adopted scheme which 
replaces the assumption of a constant M™ m . To start 
with, it is necessary to consider some indications about 
the dredge-up provided by detailed calculations of ther- 
mal pulses. They show that the maximum inward pen- 
etration of the envelope convection, over the pulse cycle, 
occurs when the surface luminosity attains its characteris- 
tic post-flash peak, L P (Wood 1981; Wood & Zarro 1981; 
Boothroyd & Sackmann 1988d). Dredge- up takes place if 
the envelope convection zone is able to penetrate into the 
carbon enriched region previously occupied by the inter- 
shell convection zone during the flash. 

In Sects. 3.1 and 3.2 below we recall and discuss the 



features which are relevant to our analysis on the occur- 
rence of the dredge-up, namely the behaviour of the peak 
luminosity Lp , and the temperature at the base of the con- 
vective envelope during the dredge-up events. In Sect. 3.3 
the method to infer M" lm for each AGB model is de- 
scribed, and the results for several values of the stellar 
mass and metallicity are presented. Section 3^ details the 
adopted inter-shell composition. 

3.1. The M c — Lp relation 

Wood & Zarro (1981) and Wood (1981) pointed out that 
at the surface luminosity maximum the outward-flowing 
luminosity, L ri is almost constant (to a fraction of a per- 
cent) throughout the envelope, from the top of the rem- 
nant inter-shell convection zone to the surface. This is 
because the thermal timescale in the envelope (i.e. for 
absorption/release of internal-gravitational energy during 
the hydrostatic re-adjustment following the flash) is much 
shorter (i.e. few years) than the evolutionary timescale 
for the envelope structure to be changed (i.e. decades to 
centuries). The constancy of L r may no longer be true 
when significant energy absorption can occur in the enve- 
lope, e.g. in stars with quite massive envelopes (M onv > 
2.5 M Q ), or large core masses (M c > 0.9 M Q ) (Wood & 
Zarro 1981). 

The fact that at the post-flash luminosity maximum 
the envelope is close to hydrostatic and thermal equilib- 
rium (Wood 1981), together with the existence of a radia- 
tive buffer over the burning shell, constitute the necessary 
physical conditions which secure the validity of the so- 
called zero-order-radiative-approximation (Schwarzschild 
1958; Eggleton 1967; Paczyhski 1970; Wood & Zarro 
1981). In brief, theory shows that Lp is essentially con- 
trolled by the core mass of the star independently of the 
outer envelope mass, in analogy with the existence of the 
M c — L relation during the quiescent inter-pulse periods 
for low-mass AGB stars. For more details on the basic con- 
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ditions that underlie the theory of a core mass-luminosity 
relationship, the reader is referred to Refsdal & Weigert 
(1970), Kippenhahn (1981) and Tuchman et al. (1983); 
see also Boothroyd & Sackmann (1988b) and references 
therein. 

In this study we adopt the M c — Lp relation presented 
by Wagenhuber & Groenewegen (1998), derived from ex- 
tensive calculations of complete AGB models (Wagenhu- 
ber 1996). It expresses the post- flash luminosity maximum 



as the sum of two terms, i.e. Lp 
where 



'(P.full) 



AL 



(P, first) j 



V.fuii) = [93000 + 2758 log(Z/0.02)] x (2) 
{M c - [0.503 + 0.82(Z - 0.02)]} 

and 

_ 1 q[1.77+2.76M c ,o-(34.0+32.4M c ,o)AM c ] /o\ 

The former term [Eq. (||)] gives the post-flash luminos- 
ity maximum (in Lq) for full- amplitude thermal pulses 
as a function of M c (in M ), with a dependence on the 
metallicity Z (in mass fraction). The latter term [Eq. (B)] 
represents a negative correction which has to be added to 
L(p f u n) in order to reproduce the fainter maximum in the 
light-curve (at given M c ) for the first pulses, expressed as 
a function of the core mass at the first He-shell flash, M Cj o, 
and the current core mass increment, AM C = M c — M c q. 



3.2. A characteristic temperature for dredge-up 

Following Wood (1981), at the post-flash luminosity peak 
the envelope convection reaches its maximum inward pen- 
etration (in mass fraction) and base temperature, T™ ax . 
At the same time, the nuclearly processed material in- 
volved in the He-shell flash is pushed out and cooled down 
to its minimum temperature over the flash-cycle, T™ m . 

To this respect, a very interesting and useful indication 
is given by Wood (1981), who reported that at the stage 
of post-flash luminosity maximum, the temperature T™ n 
always approaches a limiting constant value (logT™ n = 
6.7 ±0.1), for all core masses studied (0.57 M Q < M c < 
0.88 M Q ), regardless of the metallicity (Z = 0.001, Z = 
0.01, and Z = 0.02). 

A similar remark can be found in the work by Booth- 
royd & Sackmann (1988d) about the production of low- 
mass carbon stars. In the final discussion the authors 
stated that "it is probably no coincidence" that in all 
cases considered (i.e. full AGB calculations for stars with 
total mass and metallicity [M- t = 0.8 M Q ,Z = 0.001], 
[Mi = 1.2 M ,Z = 0.0022], [M { = 1.7 M Ql Z = 0.001], 
and [M = 2.0 M ,Z = 0.001]) the first episode of con- 
vective dredge-up is characterised by a temperature at the 
base of the convective envelope of approximately the same 
value, logT™ ax ~ 6.5. In other words, T b max - T™ in at the 
first occurrence of the third dredge-up. 

In conclusion, detailed analyses of thermal pulses sug- 
gest that the third dredge-up takes place when the base 



temperature attains or exceeds some critical value. Here- 
inafter, we will refer to it as T^ rcd . Therefore, we adopt 
the following criterion: 



no dredge-up can occur 
dredge-up occurs 



On the basis of the above considerations, it is natural to 
assume that the fulfilment of the condition 



Timax rpdrcd 
b — J b 



(4) 



for the first time during the evolution marks the onset of 
dredge-up. 

3.3. M™ ln from envelope integrations 

Thanks to the fact that at the time of the maximum pen- 
etration of the external convection the en velo pe is in hy- 
drostatic and thermal equilibrium (Sect. |3.l|) , and that 
T max > T drcd when dredge-up occurs (Sect.^J), it follows 
that the question of the minimum core mass for dredge- 
up can be analysed with the aid of a stationary envelope 
model. 

Envelope integrations are performed on the basis 
of the following scheme. According to the analytical 
demonstration of Tuchman ct al. (1983) it turns out 
that in the zero-order-radiative-approximation (see also 
Schwarzschild 1958; Eggleton 1967; Paczyhski 1970) the 
product of the temperature T c and radius R c just above 
the burning shell, depends on the surface luminosity, core 
mass M c , and chemical composition of the radiative tran- 
sition region: 



Tc_ 
10 7 



Rc 



2.291 



(5X + 3- Z) 



M c - 0.156(1 + X 



L 

To 1 



(5) 



where T c is given in K, L c , M c and R c are expressed in 
solar units, X and Z are the abundances (in mass fraction) 
of hydrogen and metals, respectively. It is worth specifying 
that the term (5X + 3 — Z) comes from the definition of 
the mean molecular weight under the assumption of a fully 
ionised gas, [i = 4/(5X + 3 — Z), and the term (1 + X) 
derives from the expression of the opacity given by the 
Thomson scattering of electrons, k — 0.2(1 + X). 

Equation (|5|) offers the boundary condition adopted 
here to determine the envelope structure of an AGB star 
at the luminosity-peak, i.e. the radius R c of the spherical 
layer below which the mass M c is contained, must coincide 
with: 



R c 



10 7 



2.291 

(3 - Z sh ) 



AL - 0.156 



10 4 



(G) 



Equation (6J) is derived from Eq. 
by Eqs. 



with L = Lp as given 
and (|3j), setting X = since in this particu- 
lar case we deal with material that has already undergone 
complete H-burning. The metallicity Z sn is that of the 
inter-shell, i.e. Z sh = X( 12 C) + X( w O) + Z ~ 0.24 ac- 
cording to the prescription given in Sect. pA. Moreover, we 



6 



Marigo, Girardi & Bressan 



assume a constant value of the temperature, T c = 10 8 K, 
which gives a reasonable order of magnitude of the min- 
imum temperature for the ignition of helium. The sensi- 



0.04 



tivity of the results to T c is discussed in Sect. 3.3.1 below 



Then, for given stellar mass, core mass, surface chem- 
ical composition, and peak- luminosity Lp, envelope in- 
tegrations are performed until we find the value of the 
effective temperature, T e g, such that M(R C ) = M c . At 
this point, the structure of the envelope is entirely and 
uniquely determined. 

3.3.1. Consistency checks 

The consistency of our prescriptions must be checked. 
First, the assumption of a constant T c = 10 8 K seems 
to be arbitrary. However, we checked that varying this 
temperature by as much as 1 dex - far more than allowed 
by reasonable changes in the He-burning reaction rates - 
produces quite negligible changes on the results of enve- 
lope integrations (e.g. the values of T™ ax , T e ff, M™ m ; see 
Table |l|). Such invariance is just expected, given the phys- 
ical decoupling between the core and the envelope as the 
M c — Lp relation holds. 

Actually, from the above discussion on the parameter 
T c , it follows that the envelope structure is also quite in- 
sensitive to the corresponding value of R c . Nevertheless, it 
is a useful exercise estimating R c in an independent way 
and comparing it with the predictions from Eq. (^|). 

Following the formalism of Tuchman et al. (1983), the 
core of an AGB star behaves essentially as a warm dwarf 
of radius: 



Ry, 



1 (M C ,T) = k(M c ,T)R cold {M c ) 



(7) 



The term warm dwarf here indicates an isothermal 
electron-degenerate stellar core with non-zero tempera- 
ture. The radius, i? C oid, of a cold dwarf (defined at zero 
temperature) is a function of its mass (and chemical 
composition) only, whereas the radius of a warm dwarf 
is larger and depends on the temperature too. The ra- 
tio i?warm/-Rcoid is expressed via the factor k{M c ,T) in 
Eq. (0), which increases with T and decreases with M c . 
However, Tuchman et al. (1983) pointed out that k typ- 
ically ranges from 1 to 3 for relevant burning shells and 
M c k, 0.6 M Q . 

The radius i? co id is calculated adopting the linear fit to 
the Chandrasekhar (1939) results, as given by Tuchman 
et al. (1983): 



>ld 



0.019 (1-0.58M C 



(8) 



where the quantities are expressed in solar units. 

Figure |l| shows R c as function of M c , calculated accord- 
ing to Eq. (p|) with T c = 10 8 K for three values of the en- 
velope metallicity, Z = 0.004, Z = 0.008, and Z = 0.019. 
The radius i? W arm as derived from Eq. (^) for four different 
values of k is also plotted. For such a choice of the param- 
eter T c , the radius R c turns out to be roughly 3 x i? co id, in 



0.03 



o.oa 



o.oi 




0.4 0.6 0.8 1 

M c (M ) 



1.2 



Fig. 1. The adopted core radius - core mass relationship 
for three different values of the metallicity Z , derived from 
Eq. (p|) with T c = 10 8 K. The radius i? W arm calculated with 
Eq. (pn is shown for few values of k, in the range 1 — 3. 
The data corresponds to the empirical white dwarf mass- 
radius relation (Koester 1987; Schmidt 1995; Koester &: 
Reimers 1996). 



agreement with the expected behaviour of warm dwarfs, 
and it is mostly comprised within 0.04 R Q and 0.02 R Q 
for the relevant range of M c in AGB stars. Finally, for 
the sake of comparison, the empirical data on the white 
dwarf mass-radius relation (Koester 1987; Schmidt 1995; 
Koester & Reimers 1996) is also shown. It is worth notic- 
ing that the observed points are mostly consistent with 
1.5 5s k < 1. This confirms the theoretical expectation 
that white dwarfs approach the zero-temperature config- 
uration of complete degeneracy. 

3.3.2. The choice of T^ red 

On the basis of the scheme just outlined, the minimum 
core mass for dredge-up of a model star of given total mass 
and surface chemical composition is found (if it exists) 
when the condition of Eq. (^j is satisfied for a unique pair 
of (Lp,T e ff) values. 

The method represents a considerable improvement 
since we can obtain useful indications about the depen- 
dence of M™ m on the envelope mass, metallicity, and peak 
strength, an important aspect so far ignored in previous 
synthetic AGB calculations (see Groenewegen & de Jong 
1993; Marigo et al. 1996a). 

There is one more prescription to be discussed, i.e. 
the adopted value of the characteristic temperature for 
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Table 2. Results from envelope integrations for a (2.0 M Q , Z = 0.004) model at the stage of the post-flash luminosity 
maximum for given (Lp,M c ) pair. Relevant quantities are shown as a function of different choices of the input 
parameter T c [see Eq. (||)]. The minimum core mass for dredge-up, M™ m , is derived assuming logT' drod = 6.4. 



M C /M Q log(L P /L ) 


log(T c /10 7 K) 


log(i?c/cm) 


log(T cff /K) 


iog(rr x /K) 


M c min /M Q 


0.5400 3.6598 


0.8 


9.6328 


3.580637 


6.48485 


0.52551 




1.0 


9.4328 


3.580642 


6.48461 


0.52556 




1.2 


9.2328 


3.580648 


6.48440 


0.52563 




1.4 


9.0328 


3.580654 


6.48408 


0.52569 




1.6 


8.8328 


3.580660 


6.48379 


0.52576 




1.8 


8.6328 


3.580667 


6.48346 


0.52584 



convective dredge-up, T 1 dred . As already mentioned, in- 
dications come from detailed analyses of thermal pulses, 
such as \ogT£ lcd ~ 6.7 in the work by Wood (1981), and 
logT drod ^ 6 5 found by Boothroyd & Sackmann (1988d). 

Wood (1981) first performed a similar kind of anal- 
ysis adopted here to infer M™ ln from envelope integra- 
tions. The author aimed at investigating the theoretical 
incapability of reproducing the observed carbon stars of 
low luminosity, i.e. no low-mass (M < 3 M Q ) and faint 
(Mbol > —5) carbon stars were predicted to form ac- 
cording to full AGB calculations. He suggested that the 
problem could be partially alleviated if the onset of the 
dredge-up episodes occurs much earlier during the TP- 
AGB evolution of low-mass stars. This implies a smaller 
M c min than found in full calculations (see Table 0). The 
question was considered by arbitrarily changing the value 
of the mixing-length parameter, initially set to a ~ 1. The 
net result was to lower systematically the minimum total 
mass for dredge-up at increasing a, whereas the corre- 
sponding Af™ m was slightly affected. Moreover, it became 
clear that, for given stellar mass M, M™ 111 increases with 
Z and, for given Z, it decreases with M. 

Boothroyd & Sackmann (1988d) also addressed the 
problem of the formation of low-mass carbon stars by 
analysing the effect of adopting high values of the mix- 
ing length parameter [e.g. a ~ 3 was necessary to produce 
a 0.81 M Ql Z = 0.002 carbon star with M c ~ 0.566 M Q 
and a luminosity at the post-flash dip log(L/L Q ) ~ 3.34, 
or equivalently M bol ~ —3.59; see also Table |j. The au- 
thors remarked that the value of a has a large effect on 
the depth in temperature of the convective envelope, but it 
has a relatively small effect on the depth in mass. In fact, 
the base of the convective envelope is located in a zone 
containing a very small mass [~ (10~ 4 — 10~ 3 )M©], but 
across which the temperature gradient is extremely steep. 
Moreover, Boothroyd & Sackmann (1988d) concluded that 
the experiment of varying a, though useful, cannot defini- 
tively solve the problem of dredge-up. 

These considerations suggest us the possibility of per- 
forming the analysis on the occurrence of third dredge-up 
from a different perspective, in the sense that we keep the 
mixing-length parameter fixed at its value (a = 1.68) de- 
rived from the calibration of the solar model, calculated 



with the complete evolutionary code of Padua (Girardi 
et al. 1996). We specify that the adopted opacities are 
taken from Iglesias & Rogers (1996; OPAL) in the high- 
temperature domain (T > 10 4 K), and from Alexander & 
Feguson (1994) in the low-tcmperature domain (T < 10 4 
K). We analyze, instead, the effect of assuming differ- 
ent values of the characteristic temperature for dredge-up. 
The lower T$ led , the smaller is the value M c min . To this 
respect, we can safely state that increasing a produces the 
same effect on M™ ln as decreasing X| drod . 

The reasons why we opt for such a choice are the fol- 
lowing. First, we want to secure the continuity of the TP- 
AGB calculations with the previous evolution, entirely fol- 
lowed using a = 1.68. This does not mean at all that 
in reality the efficiency of the convective mixing cannot 
change during the evolution of the star. However, since the 
present theoretical treatment of the convection is indeed 
inadequate to predict possible variations of a, we prefer 
to assume a conservative position. Second, the occurrence 
of convective dredge-up is not only conditioned by the 
properties of external convection (e.g. by the choice of a), 
but it also depends on the physical characteristics of the 
inter-shell region directly involved in the flash. Detailed 
calculations clearly indicate that a higher strength of the 
pulse, (usually measured by the maximum luminosity of 
the He-shell) corresponds to a greater amount of energy 
spent to expand outward the nuclearly processed material, 
that would be cooled down more efficiently (Boothroyd & 
Sackmann 1988d; Straniero et al. 1997). 

Therefore, artificially decreasing T drod is equivalent to 
hypothesising stronger thermal pulses. This is a viable ex- 
periment in a theoretical approach making use of a static 
envelope model. In fact, our analysis is not directed to 
investigate the effective penetration of the convective en- 
velope, since that would require us to carefully follow (i.e. 
with a complete stellar code) the onset and the time de- 
velopment of each thermal pulse. Obviously, an impossible 
task using only envelope integrations. The essential goal 
of this study is to simply derive useful indications on the 
possible occurrence of convective dredge- up, on the basis 
of physical aspects which can be properly handled by our 
investigative tool. 
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Fig. 2. Minimum core mass for dredge-up, M™ m , as a 
function of the current stellar mass for models with [Y = 
0.25, Z = 0.008], corresponding to three values of T drcd . 
The short-dashed line refers to the core mass at the first 
thermal pulse. 



3.3.3. The dependence of M c min on M and Z 

In this section we aim at investigating the dependence of 
M™ m on the stellar mass M and metallicity Z, for differ- 
ent values of the temperature parameter T^ rcd . 

As starting point, we explore both the cases 
logT b drcd = 6.7 (from Wood 1981) and logT b drcd = 6.5 (fol- 
lowing Boothroyd & Sackmann 1988d) to calculate M™ m 
as a function of the stellar mass for the set of models 
with original composition [Y = 0.250, Z = 0.008]. The 

results are displayed in Fig. 0, together with the case 
logT dred = g 4 

One can immediately notice that all three curves share 
a common trend, i.e. M™ m very steeply increases toward 
lower masses [M < (1.5 — 1.8)M Q ], and flattens out to 
a nearly constant value for greater masses (M > 2M ). 
The comparison of these curves with the values of the 
core mass at the starting of the thermally pulsing regime 
(short-dashed line) allows one to get an approximate es- 
timate of the core-mass increment necessary before the 
star begins to become carbon-enriched. It turns out that 
stars with masses M > 2.5M enter the TP-AGB phase 
with a core mass that is greater than the corresponding 
M™, for all the three values of temperature T dred under 
study, i.e. the third dredge-up should already occur at the 
first pulse (or, more likely, as soon as the full-amplitude 
regime is established). On the contrary, stars with lower 
masses (M < 2.5M ) are expected to evolve through a 



and M££ 



Fig. 3. Minimum core mass for dredge-up, M™ m , as a 
function of the current stellar mass for models with metal- 
licity Z = 0.004, Z = 0.008, and Z = 0.019, corresponding 
to the indicated value of logTj drcd . 

longer part of the TP-AGB without undergoing any chem- 
ical pollution by convective dredge-up. 

Moreover, it is worth noticing that each possible choice 
of T' drcd determines a different value of the minimum total 
mass, -M^j, for a star to experience dredge-up. In the Z = 
0.008 case, calculations yield ~ 1.5Af© adopting 

logT drcd = g 7) M min _ L2Mq adopting logT drcd = 6.5, 

1.1M adopting logT dred = 6.4. 
The dependence on metallicity is illustrated in Fig. |[ 
for a fixed value of logT drcd = 6.4. In agreement with 
the results from full computations of thermal pulses (e.g. 
Boothroyd & Sackmann 1988d), our analysis confirms that 
for given stellar mass, the onset of convective dredge-up is 
favoured at lower metal abundances in the envelope. It is 
worth noticing that the minimum mass, M™^, is affected 
in the same direction. According to our calculations, we 
get Af™' d ~ 1.2M for Z = 0.019, M™™ ~ 1.1M Q for 
Z = 0.008, and AP™ d ~ 1M for Z = 0.004. This result 
seems to be consistent with the observed data of carbon 
stars in the Magellanic Clouds (see Sect. ||): the faint end 
of the luminosity distribution of carbon stars in the SMC 
(with a mean metallicity Z ~ 0.004) extends to a fainter 
luminosity (M bo i ~ -2.5), than in the LMC (Af bo i ~ -3, 
with a mean Z ~ 0.008), thus suggesting a lower at 
lower metallicities. 

We like to stress once more that the onset of the mix- 
ing events in low-mass stars is a fundamental aspect, since 
it crucially affects the distribution of the carbon stars in 
the low-luminosity domain. As will be shown below (from 
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Fig. 4. Maximum base temperature, T™ ax , calculated at 
each thermal pulse (crosses connected by dotted lines) as 
a function of the core mass for a few evolving TP-AGB 
stellar models with original metallicity, Z = 0.008. Stellar 
masses at the onset of the TP-AGB phase are indicated 
in solar units nearby the corresponding curves. The solid 
horizontal line marks the critical temperature (log = 
6.4) for the occurrence of the third dredge-up. 



a systematic analysis), the value logX^ rcd = 6.4 is found 
suitable to fulfill the above constraint (Sect. ^) in the Mag- 
ellanic Clouds. 

It is worth specifying that the results illustrated in 
Fig. U give M™ m as a function of the current total mass 
(not the mass at the first thermal pulse). This is partic- 
ularly relevant in the low-mass domain, where the curves 
become very steep, i.e. minor variations of the total mass 
correspond to quite different values of Af™ m . Since TP- 
AGB calculations are performed with the inclusion of mass 
loss by stellar winds, it follows that in the case of low-mass 
stars a small reduction of the envelope mass may signifi- 
cantly delay, or even prevent the onset of dredge-up. 

Figure ^ shows the maximum temperature at the base 
of the convective envelope, T™ ax , at the stage of post-flash 
luminosity peak, as a function of the core mass during 
the evolution of a few stellar models with different masses 
M at the beginning of the TP-AGB phase, and original 
metallicity Z 



0.008 (see Sect. 3.5 for further details 



of the models). The adopted temperature parameter is 
ydred = 64 

Envelope integrations indicate that if the envelope 
mass is too small, the star fails to attain the minimum 
temperature T^ rcd , whatever the core mass is. See, for in- 
stance, the track for the M = 0.8M Q star. The model with 



M = 1.2M corresponds to the lowest mass for a star to 
become carbon-enriched. Higher mass stars are able to at- 
tain hot base temperatures, so that they are expected to 
experience dredge-up already since the first pulses. 

The present method also provides the possibility of 
predicting the stage beyond which dredge-up cannot oc- 
cur any longer during the final TP-AGB evolution. The 
shut-down of dredge-up may happen when the envelope 
mass becomes so small that the condition 7 1 ™ ax > ydred - 1S 
no longer satisfied. This circumstance occurs for the last 
pulses of the tracks displayed in Fig. ||, characterized by 
the drastic reduction of the envelope mass due to stellar 
winds. This point represents another important improve- 
ment for synthetic TP-AGB models, since the end of the 
dredge-up episodes is not fixed a priori (a necessary as- 
sumption if one adopts a constant value for M™ ln ), but 
naturally derives from calculations. 

3.4- The inter-shell chemical composition 

The remaining ingredient of our models is the composi- 
tion of the dredged-up material. We use the same formal- 
ism as in Marigo et al. (1996a), to whom we refer for an 
extensive description of this point. Suffice it to recall that 
our prescription is based on the detailed calculations of 
the inter-shell nucleosynthesis carried out by Boothroyd 
& Sackmann (1988c). They find that the convective inter- 
shell after a thermal pulse and just before the penetra- 
tion of the envelope is composed by 20 — 25 % of 12 C, 
~ 2 % of 16 O, and the remainder of 4 He (in mass frac- 
tions). These proportions are reached after the first few 
thermal pulses, and do not significantly depend on the 
core mass and metallicity. 

Recently, Herwig et al. (1997) have pointed out that 
the abundances in the inter-shell will differ significantly 
from the results of standard calculations (i.e. without 
extra-mixing as in Boothroyd & Sackmann 1988c), if in- 
stead one applies their overshooting scheme to all con- 
vective zones. The resulting chemical distribution of the 
dredged-up material typically consists of 50% in 12 C, 25% 
in 16 0, and 25% in 4 He. 

In this study we opt for the conservative standard pre- 
scription of Boothroyd & Sackmann (1988c). The follow- 
ing fixed abundances are adopted: 22% for 12 C produced 
by the triple alpha reaction and 2% for 16 produced by 
the reaction 12 C(a, 7) 16 0. 

3.5. Other computational details 

The evolutionary calculations presented in this study 
have been carried out for a dense grid of stellar masses 
(0.8M Q <j M\ < 5Mq) and two values of the initial metal- 
licity, Z = 0.008 and Z = 0.004. In addition to the new 
treatment of the third dredge-up already presented, the 
present calculations are based on a synthetic TP-AGB 
model, having the same general structure as described 
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in Marigo et al. (1996a), with some major revisions and 
updates illustrated in Marigo et al. (1998) and Marigo 
(1998ab). The reader should refer to these works for all the 
details. Suffice it to recall here that the synthetic model 
combines analytical relations with a complete stellar enve- 
lope, which is integrated from the photosphere down to the 
core. The evolution of a TP-AGB star is followed from the 
first thermal pulse up to the complete ejection of the enve- 
lope, adopting the prescription for mass loss as suggested 
by Vassiliadis & Wood (1993). A homogeneous set of ac- 
curate analytical relations (e.g. the core mass-luminosity 
relation, the core mass-interpulse period) is taken from 
the recent work by Wagenhuber & Groenewegen (1998). 
Finally, it is worth remarking that the occurrence of en- 
velope burning (or hot-bottom burning) in more massive 
TP-AGB stars (M ^ 3.5M ) is consistently followed to 
account for the possible break-down of the M c — L relation 
(see Marigo 1998b). 



4. The observed carbon star luminosity functions 
in the MCs 

The Magellanic Clouds - LMC and SMC - are an ideal 
site for testing evolutionary theories of AGB stars, thanks 
to the fortunate combination of several factors: 

— The extinction to both Clouds is low and relatively 
uniform. 

— The galaxies are far enough (~ 50 kpc to the LMC, 
~ 63 kpc to the SMC) that their depth along the line 
of sight can be neglected to first order, and the same 
distance can be assigned to all the stars. 

— At the same time the galaxies are near enough to al- 
low the resolution and analysis of single stellar sources 
by means of both photometric and spectroscopic tech- 
niques. 

— Both Clouds have a considerable number of populous 
star clusters whose age (Elson & Fall 1985; Girardi 
et al. 1995) and chemical composition (Olszewski et 
al. 1991) can be estimated with a good accuracy, and 
which contain a considerable number of AGB stars 
(Frogel et al. 1990; Westerlund et al. 1991). 

— A large amount of observed data on AGB stars is nowa- 
days available (see the reviews by Westerlund 1990 and 
Groenewegen 1997), thanks to large scale surveys rang- 
ing from the optical to the radio domain, and spectro- 
scopic analyses for the estimation of stellar parameters 
(e.g. effective temperature, mass-loss rate, wind veloc- 
ity, surface gravity, surface chemical abundances). 

For the observed carbon star luminosity function in the 
LMC we make use of the recent compilation of Costa & 
Frogel (1996). The data refer to 895 carbon stars, selected 
among the original sample of 1035 stars found by Blanco 
et al. (1980) and Blanco & McCarty (1983) by means of 
low-dispersion red grism spectroscopy in 52 fields of the 
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Fig. 5. The observed LFs of carbon stars detected in the 
fields of the LMC (solid line) and SMC (dashed line), nor- 
malised to the total number of stars in the samples. The 
data for the LMC are taken from Costa & Frogel (1996), 
adopting a distance modulus fj-o = 18.5. The data for the 
SMC is as derived by Groenewegen (1997) from the data 
compiled by Rebeirot et al. (1993), with a distance mod- 
ulus no = 19.0. 



galaxy, distributed over an area of 8° x 8°. The identi- 
fication of carbon stars was based on the recognition in 
stellar spectra of a typical blend of the CN bands at near- 
infrared wavelengths (AA7910, 8100, 8320). It is worth re- 
marking that this sample of carbon stars can be considered 
almost complete and not limited by apparent magnitude, 
since the faintest carbon stars are at least 2 mag brighter 
(I ~ 15) than the faint limit of detection (I ~ 17). More- 
over, the possible contribution of high luminosity carbon 
stars, which would have been missed in the survey because 
of dust obscuration, is likely to be rather small (^ 3% for 
Mboi < —6 according to Groenewegen & de Jong 1993). 
Costa & Frogel (1996) performed RI and/or JHL pho- 
tometry of the selected 895 carbon stars. The apparent 
bolometric magnitudes raboi were estimated first applying 
proper bolometric corrections to the K magnitude of a 
subsample of stars with available JHK photometry. On 
the basis of these results the authors derived a linear re- 
lationship, expressing (Iq — mbol) a s a function of the red 
colour (R — J)o, which was then used to obtain mbol for 
the majority of carbon stars with only RI data. The au- 
thors quoted a mean error on the evaluation of Wboi of 
approximately ±0.34 mag. 

The observed CSLF in the LMC is plotted in Fig. | 
(solid line). The absolute bolometric magnitudes, Mbol) 
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are obtained from the apparent ones adopting a distance 
modulus = TO boi — Afboi = 18.5, as indicated by West- 
erlund (1990) and favoured by several recent works (e.g. 
Wood et al. 1997; Oudmaijer et al. 1998; Madore & Freed- 
man 1998; Salaris & Cassisi 1998; Panagia 1998). 

The resulting CSLF is quite similar to those obtained 
from previous samples of carbon stars in the fields of the 
LMC (Westerlund et al. 1978; Richer et al. 1979; Cohen 
et al. 1981). It extends approximately from Afboi = —3 up 
to Mboi = —6.5, with the peak located at around Mboi = 
—4.875 (centre of the bin). The dispersion of the observed 
distribution is estimated to be about a = 0.55 mag. 

From a comparison with the CSLF in selected clusters 
of the LMC (Frogel et al. 1990) it turns out that the two 
distributions are quite similar, except for few subtle differ- 
ences, namely, the presence of somewhat more extended 
faint and bright tails in the field distribution relative to 
that of the clusters. According to Marigo et al. (1996b), 
the explanation of this resides in the presence of two dif- 
ferent periods of quiescence in the cluster formation his- 
tory of the LMC, shaping the age (and progenitor mass) 
distribution of C stars. These gaps are found at ages of 
- (3 - 12) x 10 9 yr, and - (2 - 6) x 10 s yr (see also Gi- 
rardi et al. 1995), which roughly correspond to the ages of 
the least (~ 1.2M©) and most (~ 4M Q ) massive carbon 
star progenitors, respectively. 

Figure || also shows the observed CSLF in the SMC 
(dashed line), as derived by Groenewegen (1997) from a 
sample of 1636 stars observed by Rebeirot et al. (1993). 
The adopted bolometric correction are those by Wester- 
lund et al. (1986). The adopted distance modulus for the 
SMC is hq = 19.0 (see Westerlund 1990 for a review of 
SMC distance determinations). 

The first bin of this CSLF contains all stars fainter 
than Afboi = —2.5. They are less than 3 % of the total 
sample of SMC carbon stars. Their intrinsic magnitudes 
can be as low as Mboi = —1.8, as found by Westerlund 
et al. (1995). The formation of carbon stars at such low 
luminosities can hardly be understood as the result of the 
third dredge-up in single AGB stars. Most probably, they 
arise from the evolution in close binary systems (see e.g. de 
Kool & Green 1995; Frantsman 1997), in which the surface 
of a low-mass star has been contaminated by the ejecta of 
a former AGB star companion. The low-luminosity car- 
bon stars will not be considered in this work, since we 
are dealing with single star evolution. For all practical 
comparisons, we will consider the CSLF in the SMC as 
extending down to Mboi = —2.5. 

The CSLFs are significantly different in the two Magel- 
lanic Clouds. In the SMC the CSLF is much broader, and 
presents a fainter peak location (at Afboi ~ —4.375; centre 

of the bin) than that in the LMC (Af bo i 4.875). This 

indicates that LMC carbon stars are, on the mean, about 
60 % brighter than their SMC counterparts. As remarked 
by Groenewegen (1997), explaining this difference quan- 



titatively represents a challenge for synthetic TP-AGB 
models. 

5. The theoretical carbon star luminosity function 

In order to simulate the observed CSLF in the Magellanic 
Clouds, we employ the general scheme already used in 
Marigo et al. (1996a). For the sake of completeness, we 
here simply recall the basic steps. 

Dividing the luminosity range of interest (e.g. —2.5 > 
Afboi > —7.0) into a suitable number of intervals, for each 
k th luminosity bin of width A Afboi, we are required to eval- 
uate the number Nk of carbon stars which are expected 
to be in transit at the present epoch: 

M max 

N k = / A^AfOdAfi 



M. ml 



$(Af i ,t)Ar fc (Af i )dAf i 



(9) 



where $(Mj,f) is the stellar birth rate, i.e. the number of 
stars formed with mass M\ at the epoch t, per unit of mass 



and time (see Sect. 5.1); and ATk(M{) is the time each star 
of mass Af; spends in the k th luminosity interval. Thus, 
Nk(Mi) expresses the number of carbon stars evolved from 
a progenitor with initial mass Afj and a current (i.e. at the 
present time) luminosity inside the k th bin. 

The integration is performed over the mass range 
[Afi min , Afi max ]. We denote with Afi min the maximum of 
the lower mass limit for the formation of carbon stars and 
the mass of the star whose lifetime is just the age of the 
galaxy (~ 0.87Af Q for T G = 15 Gyr). We refer to Af; max as 
the minimum of the upper mass limit for carbon star for- 
mation, and the maximum mass of stars which experience 
the AGB phase, Af up . In our case, the evolution of the 
progenitor stars has been calculated with the adoption of 
an overshoot scheme applied to the convective zones, from 
the ZAMS up to the starting of the TP-AGB (Girardi et 
al. 1998). The overshooting parameter for core convection 
adopted for intermediate-mass stars (A c = 0.5) charac- 
terises a moderate efficiency for this process. This fixes 
Af up ~ 5Af Q . Classical models without overshoot would 
predict a higher value, e.g. Af up ~ (7 — 8)M Q . 

Let us now discuss the meaning of the two basic func- 
tions in Eq. @, namely $(Af i; i) and Ar fc (Afi). 

5.1. The stellar birth rate 

It is useful to split the stellar birth rate <I> into the product 
of two functions (see Tinsley 1980): 

*(Mi,t)=^(Mi)^(t) (10) 

where 4>(M\) is the initial mass function (IMF) express- 
ing the fraction of stars forming with masses in the range 
[Mi, Af; + dAfi]. This latter is normalised so that: 



Af; 0(Af;) dM = 1 



(11) 
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In this work we adopt the classical power-law IMF of 
Salpeter (1955), <j)(M{) oc Mr (1+x] , with x = 1.35. The 
function ip(t) is the star formation rate (SFR) and rep- 
resents the mass of primordial gas that is converted into 
stars per unit of time at the epoch t. 

In the particular case under consideration [see Eq. (||)], 
it is required to evaluate the number of carbon stars with 
initial mass M; and lifetime t* which are born at the epoch 
t = Tq — t* (and hence are presently observable), where 
Tq is the age of the galaxy. This can be expressed by: 



$(M i; Tq - T*) = 0(Mi) V[T G - r* (Mi)] 



(12) 



Since there is a unique correspondence between the life- 
time of a star (of given metallicity) and its initial mass, 
t* = T*(Mi), it is possible to consider $ as a function 
solely of the stellar mass. 

5. 2. The rate of brightening 

The quantity Arfe(Afj) is the time that a star with ini- 
tial mass Mi spends in the k th luminosity bin of width 
AMboi- It is directly derived from the TP-AGB evolu- 
tionary tracks. 

A useful approximation for interpreting the results is 
that the quiescent luminosity of TP-AGB stars follows an 
almost constant rate of brightening (see e.g. Iben 1981): 



dMbo 
dr 



constant 



(13) 



Actually, this relation holds if the following conditions are 
met. First, the relation between the quiescent luminosity L 
and the core mass M c should be linear. This would corre- 
spond to consider only the terms (5a) and (5b) of Eq. (5) 
giving the M c — L relation adopted by Marigo (1998b). 
This implies the neglect of both the initial sub-luminous 
evolution of the first inter-pulse periods, and the possible 
deviations towards higher luminosities caused by envelope 
burning. This assumptions are reasonable for carbon stars 
since their formation generally occurs after the onset of the 
full amplitude regime, and only when envelope burning is 
not efficiently operating. In this case: 



dM c 
dL 



(14) 



On the other hand, the evolutionary rate (i.e. the growth 
rate of the core mass) is proportional to the stellar lumi- 
nosity: 



dM c v L 



(15) 



In the above equations K\ and are suitable constants. 
If the changes in the envelope abundance of hydrogen, X, 
are small enough, it follows that dM c can be eliminated 
from the above equations giving dlnL/dr — constant, 
which is equivalent to Eq. (|l^) . This equation tells us that 



the rate of brightening along the TP-AGB is nearly con- 
stant, regardless of the core mass, the total mass, and the 
luminosity level. 

In our computations, the rate of brightening for carbon 
stars is found to be constant within ~ 35 % at most. Typi- 
cal values for this quantity are ~ —(3.5 — 7) 10~ 7 magyr _1 . 
Under these conditions, it follows that nearly the same 
number of stars is expected to populate each magnitude 
interval along the AGB phase of coeval (e.g. a cluster) 
stellar populations. If we consider only the evolution of 
the quiescent luminosity (see, however, Sect. 5A below), 
the resulting CSLF would be described by a box-like dis- 
tribution. 

A final remark is worth making. Since dredge-up peri- 
odically reduces the core mass, then it is expected to affect 
the average evolutionary rate dM c /dr. This implies that 
Eq. ( |l5| ) is valid only if either dredge-up is not operating, 
or A = constant as in our prescriptions. In contrast, for A 
variable with time, neither K2 nor the rate of brightening 
dMboi/dr would be constant. 

5.3. Transition luminosities 

The calculation of the CSLF requires also the knowledge 
of: 

— the luminosities marking the transition from the M- 
class (C/O < 1) to the C-class (C/O > 1), 

— the maximum luminosities at which the star is still a 
carbon star, 

as a function of the initial mass. 

Figures ^| and [7] show the transition luminosities as 
a function of the initial mass for the Z = 0.008 and 
Z = 0.004 cases, respectively. The choice of dredge-up 
parameters adopted to construct these plots will be dis- 
cussed later in Sects. || and 0. Several interesting features 
can be noticed from the inspection of these figures. 

First, the formation of carbon stars may be prevented 
because of missing or insufficient dredge-up in the range 
of lower masses, or significantly delayed due to the oc- 
currence of envelope burning in the most massive AGB 
stars. 

Second, for the majority of low-mass carbon stars 
(with no envelope burning; Mi < 3.5M Q ; see Marigo 
1998b) the maximum luminosity just coincides with the 
tip of the AGB reached at the end of their evolution. On 
the contrary, in more massive carbon stars (Mi > 3.5M Q ) 
the occurrence of envelope burning can lead to a different 
configuration. It may happen that the onset of envelope 
burning takes place after the star has already experienced 
a sufficient number of dredge- up episodes to become a car- 
bon star. Nuclear burning at the base of the envelope sub- 
sequently destroys carbon, lowering the C/O ratio again 
below unity. Then, the star is expected to evolve like an 
oxygen- rich star, until a new transition to the C-class pos- 
sibly occurs in the very late stages of the evolution. In fact, 
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Fig. 6. Transition bolometric magnitudes vs. initial 
masses of TP-AGB stars with original metallicity Z = 
0.008. Regions labelled with letters A, B and C refer to 
three choices of the dredge-up parameters. See the text for 
more details. 

after envelope burning has extinguished, the star eventu- 
ally undergoes further convective dredge-up of carbon be- 
fore leaving the AGB phase (see Frost et al. 1998). Then 
the star would appear as a carbon star in two distinct evo- 
lutionary stages, being oxygen-rich in between. This latter 
configuration can be clearly noticed in Figs. ^ and [?]: at 
masses higher than about 3.5M Q , a thin strip indicates 
the domain of the high-luminosity carbon stars, i.e. those 
which experienced a late conversion to this spectral type. 

In this case, the transition to the carbon star phase 
may occur under conditions of severe mass-loss. Thus, 
the brightest carbon stars are expected to be obscured 
by thick layers of ejected material, then not detectable 
at visual wavelengthes. Actually, high-luminosity obscured 
carbon stars have been recently discovered in the LMC, 
by means of infra-red observations (van Loon et al. 1998). 
Since they are not included in the observed CSLFs in the 
Magellanic Clouds (Sect. |J), we will not consider them 
here. The formation of these carbon stars will be inves- 
tigated by means of the present models in a forthcoming 
paper. 

5.4- The effect of the low-luminosity dip 

The fact that low-mass (M; < 3M ) TP-AGB stars un- 
dergo a long-lived subluminous stage, soon after the oc- 
currence of a thermal pulse, plays a crucial role in de- 
termining the extension of the low-luminosity tail of the 
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Fig. 7. Transition bolometric magnitudes vs. initial 
masses of TP-AGB stars with original metallicity Z = 
0.004. See the text for more details. 

CSLF (Boothroyd & Sackmann 1988ad). It turns out that 
without accounting for this effect the theoretical distribu- 
tion cannot extend down to the observed faint end (see 
also Groenewegen & de Jong 1993), unless quite extreme 
(and unlikely) assumptions are adopted, e.g. the occur- 
rence of very efficient (A ~ 1) dredge-up events since the 
first weak thermal pulses. On the contrary, the effect of 
the post-flash luminosity maximum on the CSLF can be 
neglected, due to its very short duration (Boothroyd & 
Sackmann 1988a). 

In our calculations, the drop in luminosity below the 
M c — L relation is included in a very simple way. On 
the basis of the results from complete calculations of 
thermal pulses (Boothroyd & Sackmann 1988a) we con- 
struct a function, P(L — AL), expressing the probability 
that a star with quiescent luminosity L, predicted by the 
M c — L relation, may be found at a lower luminosity level, 
L — AL. The maximum excursion towards fainter lumi- 
nosities is assumed to be AlogL = 0.5, corresponding to 
AMbo\ = 1-25 mag. The probability function is then esti- 
mated for a suitable number of luminosity intervals (5 in 
our case), with a width equal to the adopted bin resolu- 
tion, i.e. Alogi = 0.1, or equivalently AM DO i = 0.25. 

The probability distribution, normalized to unity, is 
assumed as following: 

Fo(logi) = 0.35 

Pi(logL-O.l) = 0.28 

P 2 (logL-0.2) = 0.20 

P 3 (log£-0.3) = 0.10 (16) 
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P 4 (logL-0.4) 
P 5 (logL-0.5) 



0.05 
0.02 



Then, the number Nk of carbon stars [see Eq. (JsJ) ] pre- 
dicted to populate the k th luminosity bin is calculated 
according to: 



N k = N k n P + 



E 

i=l 



(17) 



where Nk,o Po is the number of objects, following the 
M c — L relation, which are expected to be presently cross- 
ing the fc th interval. The term N^+i Pi represents the num- 
ber of stars having higher quiescent luminosity within 
the (k + i) th bin, but presently contributing to the fc th 
bin, because of the overlap with the corresponding low- 
luminosity dip. In other words, the luminosity function 
derived just applying the M c — L relation, is then con- 
volved with the function P(L) which mimics the shape of 
the sub- luminous stages driven by thermal pulses. 

We have actually adopted the same prescription for 
all the stars considered, although the duration and the 
depth of the post-flash luminosity dip are less pronounced 
in stars of higher core masses (Boothroyd & Sackmann 
1988a). Nevertheless, we expect that our approximation 
does not compromise the validity of the results. This be- 
cause most of the carbon stars are found in the mass range 
Mi < 2.5-M©, for which well extended luminosity dips, 
consistent with Eq. (|l7|), are expected. Moreover, it will 
be shown in Sect. ||that these stars completely determine 
the faint end and peak location of the theoretical CSLFs. 
More massive AGB stars, instead, are predicted to have 
less-pronounced luminosity dips. However, they populate 
mostly the bright wing of the CSLF, which is affected by 
several other factors (such as the efficiency of envelope 
burning, mass-loss, and the SFR; see Sect. ||). 

6. The calibration of the dredge-up parameters 
for the LMC 

The dredge-up parameters, Tj dred and A, need to be cal- 
ibrated on the basis of some observational constraint. 
First, we aim at reproducing the CSLF in the LMC, 
which has already been the subject of previous works 
(e.g. Groenewegen & de Jong 1993; Marigo et al. 1996a). 
To this aim, we employ the results of synthetic calcula- 
tions for TP- AGB stars with initial chemical composition 
[Y = 0.250; Z = 0.008]. The adopted metallicity is in 
agreement with the observational determinations for the 
LMC, specially with regard to the young components (e.g. 
Russell & Dopita 1990; Dopita et al. 1997). Moreover, Ol- 
szewski et al. (1991) and de Freitas Pacheco et al. (1998) 
find similar metallicitics in the LMC star clusters younger 
than about 4 Gyr. 

By means of the method outlined in Sect. ^, we have 
calculated the CSLF for different combinations of the two 




Fig. 8. Theoretical luminosity functions of carbon stars 
(solid lines) for three different combinations of the dredge- 
up parameters as indicated. The observed distribution in 
the LMC (dotted line) is plotted for comparison. 



parameters, while keeping fixed the other input prescrip- 
tions (namely the IMF and SFR). Specifically, we em- 
ployed Salpeter's law (1955) for the IMF, and the history 
of SFR suggested by Bertelli et al. (1992) in their study of 
the CMDs of field stars in selected areas of the LMC. The 
SFR has been moderate from the beginning up to about 
5 Gyr ago and since then a factor of ten stronger (see the 
top panel of Fig. [To]). The sensitivity of the results to var- 
ious possible choices of the SFR as a function of the age, 



and to the assumed IMF, are explored in Sect. 3.3. 

In Fig. H we show the predicted LF of carbon stars for 
three pairs of dredge- up parameters. The values are the 
following: 

logT^ rcd = 6.5, A = 0.65 



logT drod = g 4) A 

logT drod = g_ 4j A 



0.65 
0.50 



The corresponding transition luminosities from the 
oxygen-rich to the carbon-rich domain are plotted as a 
function of the initial mass in Fig. [| It is apparent that 
the formation of carbon stars of lower masses is favoured 
at decreasing ed and increasing A. On the contrary, the 
luminosity at the tip of the TP- AGB seems to be essen- 
tially independent of the dredge-up parameters (at least 
for the cases under consideration). A reliable calibration 
of these parameters requires one to disentangle which dis- 
tinctive effects are produced by one parameter without the 
influence of the other. This because it would be desirable 
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Table 3. 



Minimum initial mass, M^f^, of carbon stars' 



0.9 



0.8 



0.7 



0.6 



0.5 



— e — Z = 0.008, logT^ ed = 6.5 
— * — Z = 0.008, logT^ 6 " = 6.4 
- Z = 0.004, logT dred = 6.4 




3 
(M ) 



Fig. 9. Core mass at the first dredge-up event as a func- 
tion of the stellar mass at the onset of the TP-AGB phase, 
for two values of the temperature parameter J^ red , and 
metallicities Z = 0.008 and Z = 0.004. 



to derive a unique pair of proper values, or at least to re- 
duce the suitable space of parameters as much as possible. 



6.1. The temperature parameter. 



1 b 



Let us first discuss the role of T^ Ied 



We recall that this 
parameter determines the minimum core mass for dredge- 
up, M™ m , as a function of the stellar mass and metallicity. 

In a previous analysis (Marigo et al. 1996a) the repro- 
duction of the CSLF in the LMC was obtained adopting 
M" lm — O.58M as a constant parameter, for all stellar 
masses and chemical compositions (see also Groenewegen 
& de Jong 1993 for similar results). In contrast, we have 
shown in Sect. 3.2 that with the new algorithm based 
on envelope integrations, at fixed T drod the behaviour of 
M™ m as a function of M is far from being constant, but 
steeply increasing for lower values of the mass. Moreover, 
the onset of the third dredge-up is favoured at lower Z. 

Figure [)] illustrates the core mass - necessarily > Af™ 111 
- just before the occurrence of the first dredge-up episode 
as a function of the total mass at the starting of the TP- 
AGB phase for the evolutionary models calculated with 
initial metallicity Z = 0.008 and Z = 0.004. The compar- 
ison with the previous Fig. shows the following differ- 
ences. 

— First, the core mass at the first dredge-up event does 
not increase so steeply as M" lm in the range of low- 
mass stars M < 2.5M©. This is the result of the ex- 
treme sensitivity of M™ m to the current stellar mass. 



progenitors, and age of the relative stellar generation for 
different values of the temperature parameter T' dred and 
original metallicities. The corresponding mass at the onset 
of the TP-AGB is also indicated inside parentheses. 



z 


logT b drcd 


M carb 
J,J mm 


Age (yr) 


0.008 


6.7 


2.00 (2.0) 


1.211 x 10 9 




6.5 


1.83 (1.8) 


1.514 x 10 9 




6.4 


1.32 (1.2) 


3.890 x 10 9 


0.004 


6.4 


1.13 (1.1) 


5.220 x 10 9 



Even small reduction of the envelope actually prevents 
the formation of carbon stars in the steepest part of 
the curve shown in Fig. || 

— Second, the core mass at the onset of the third dredge- 
up is found to increase for stars with M PS 2.5M©. 
In this mass range, it reflects the behaviour of the 
core mass at the first thermal pulse (which increases 
with stellar mass), and not the flattening predicted for 

M min 

— Finally, the core masses in the low-mass range indi- 
cated in Fig. H are somewhat larger than the corre- 
sponding ones shown in Fig. |[ The reason is partly 
due to the fact that the results illustrated in this lat- 
ter plot are derived from calculations in which Lp = 
L( P) f u n) + AL(p i g rst ) (see Sect. 3.1), whereas the term 
Ai(p first ) has been neglected in the envelope integra- 
tions performed to construct Fig. ||. This simplification 
was adopted because AL(p i g rs t) contains a dependence 
on the core mass at the first thermal pulse [see Eq. (||)] , 
whereas at that stage we simply aimed at deriving the 
behaviour of M ™ m as a function of the current stellar 
mass. The neglect of Ai(p ) g rst ) leads one to overesti- 
mate Lp for a given core mass, thus yielding slightly 
lower values of M" lm . 



Comparing the curves for the same logTj drcd = 6.4 
and metallicities Z = 0.008 and Z = 0.004, it turns 
out that at decreasing metallicity the core mass at the 
first dredge-up episode is smaller in the low-mass range 
(M £3 2.2M0), whereas it is larger at higher stellar masses 
((M > 2.2M©). This is explained as following. The fact 
that at lower metallicities the minimum base tempera- 
ture for dredge-up is attained at lower values of Lp (and 
hence M c ) is only relevant in the low-mass domain, since 
the onset of the third dredge-up at larger stellar masses 
essentially reproduces the trend of the core mass at the 
first thermal pulse, which actually increases at decreasing 
metallicities for given stellar mass. 

Moreover, if we now compare the two curves plotted in 
Fig. |, for two values of T drod and and the same Z = 0.008, 
it turns out that the temperature parameter determines 
the minimum mass for a TP-AGB star to experience chem- 
ical pollution by convective dredge-up. The lower is T dred 
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the lower is this limit mass. Then, for a given mass-loss 
prescription applied to the previous evolution (i.e. during 
the RGB and part of the AGB up to the onset of dredge- 
up), we can derive the corresponding stellar mass at the 
ZAMS. In other words, T drod fixes a lower limit to the min- 
imum initial mass, M^f, for a star to become a carbon 
star, and hence an upper limit to the oldest age of the 
stellar population possibly contributing to the observed 
distribution (see Table ||) . 

Another interesting point to be stressed is that the 
results shown in Fig. |^ are determined by the temperature 
parameter only (since A cannot have any effect prior to 
dredge-up). This fact allows us to predict already that 
the value logT drcd = 6.5 should be ruled out, otherwise 
we would miss the low-luminosity tail of the LMC CSLF. 
Let us consider the minimum (with M c = Q.§77Mq) of the 
corresponding curve in Fig. ^, referring to a star with M = 
2.5M Q . Even in the most favourable (but quite unlikely) 
hypothesis, that the transition to the C-class occurs soon 
after the first dredge-up episodes (this would require A ~ 
1), the corresponding luminosity derived from the M c — L 
relation would be log(L/L Q ) ~ 3.706. This is equivalent 
to Mboi ~ —4.564, and even adding 1.25 mag (i.e . the 
typical depth of the low-luminosity dip; see Sect. 5.4), we 
get a value that is more than ~ 0.3 mag brighter than 
the observed faint end of the distribution for the LMC, at 
Mboi ~ -3. 



In the case of log 7]' 



drcd 



6.4, the minimum of the 



curve is at M = 1.8M Q , with M c = O.538Af . Applying 
the same reasoning, we would obtain the faintest carbon 
stars at Mboi ~ —2.795. The latter value is consistent with 
observed faint end of the CSLF in the LMC. 

Even if a full discussion on the carbon star transition 
luminosity as a function of the stellar mass cannot be de- 
veloped without considering the combined effect of T dred 
and A, we can conclude from the previous analysis that 
the temperature parameter, T drcd , mainly determines the 
faint end of the CSLF. 



6.2. The efficiency parameter: A 

Let us now analyse the role of the dredge-up efficiency pa- 
rameter A. We remind that A = 0.65 was the suitable value 
estimated by Marigo et al. (1996a) such that, together 
with M c min = 0.58M Q , the CSLF in the LMC was repro- 
duced. The top and middle panels of Fig. || show the the- 
oretical distributions of carbon stars adopting A = 0.65, 
for logT drcd = 6.5 and logT drod = 6.4, respectively. None 
of them fits the observed histogram. 

In the case with (logT drod = 6.5, A = 0.65), the theo- 
retical distribution is systematically shifted to brighter lu- 
minosities, with respect to the observed one. The opposite 
situation is found with the combination (logT drcd = 6.4, 
A = 0.65), which produces the peak of the predicted CSLF 
at a fainter magnitude than observed, with an excess of 
carbon stars at lower luminosity, and conversely a deficit 



in the brighter domain. Note, however, that contrary to 
the previous result with (logT drcd = 6.5, A = 0.65), both 
extremes of the observed distribution are fitted. 

This circumstance suggests that the we are only re- 
quired to shift the peak of the histogram into the right lu- 
minosity bin. We find that this can be obtained by decreas- 
ing the parameter A while keeping fixed logT drcd = 6.4, 
rather than varying the temperature parameter itself. 
There are two main reasons. First, the constraint set by 
the faint end of the distribution, which has been already 
fulfilled with the choice logT drcd = 6.4, would be violated 
with a higher value of T dlcd . Second, we expect that a 
suitable decrease of A could move the peak into the right 
location, and reproduce the faint tail of the distribution 
at the same time. In fact, the transition luminosity for 
a star with a very low-mass envelope (~ IM@) is only 
slightly affected by changes of A (provided that A <; 0.4), 
since 2 — 3 dredge-up episodes generally suffice to cause 
the transition to carbon star. 

Concluding this discussion, we remark that once T drcd 
has been determined, the efficiency parameter A essen- 
tially controls the peak location of the CSLF. We find that 
the best fit to the observed CSLF in the LMC is obtained 
adopting logT drcd = 6.4 and A = 0.5. This fit is shown in 
the bottom panel of Fig. [§. 

6.3. Sensitivity to the SFR and IMF 

The relative weight of different stellar generations con- 
tributing to the observed CSLF is determined by the SFR 
history in the host galaxy, and by the IMF. In this section, 
we investigate the sensitivity of the predictions on both 
functions. Let us first analyse the effect due to the SFR 
history, while fixing the IMF equal to Salpeter's (1955) 
prescription. Then, we explore the sensitivity of the CSLF 
to different choices of the IMF. 

We now address the following questions: 

— To what extent is the CSLF sensitive to the underlying 
SFR ? 

— Are there any distinctive features of the CSLF setting 
constraints on the history of SFR ? 

— What is more important in shaping the observed 
CSLF, the efficiency of the dredge-up process or the 
time evolution of the SFR ? 

First of all, it is clear that the possible indications the 
CSLF may give on the history of the SFR are restricted 
to the age interval defined by the minimum and maxi- 
mum initial mass of the progenitors (M^ b and M^J, 
respectively). Specifically, we deal with the age interval 
ranging from about 3.9 x 10 9 yr down to 1.1 x 10 s yr, 
corresponding to stellar lifetimes for M£5J~ = 1.32M Q 
and Mg£~ = 5M Q , respectively. The lower limit of the 

mass range is determined by the temperature parameter, 
logT ma X = 6 4 ( gee g ect Q 
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Fig. 10. Schematic representation of three different pre- 
scriptions for the SFR history in the LMC (in arbitrary 
units). In the top panel the law suggested by Bertelli et al. 
(1992) is shown [ip(t) — tpi]. See the text for more details. 



An interesting result comes out when comparing the 
theoretical distributions, each corresponding to a single 
stellar generation, as a function of the initial mass of the 
progenitor, or equivalently, of the age. These theoretical 
CSLFs are shown in Fig. [nj for the case logT^ rcd = 6.4 
and A = 0.5. For each stellar mass (and hence age), the 
CSLF has a shape determined essentially by the convolu- 
tion of a box-like function (see Sec t. |5.2| ) with the prob- 
ability function described in Sect. 5.4. From this figure, 
it turns out that all the histograms relative to stars with 
initial masses Mj < 2.5M have a faint wing extending 
within the same luminosity range, and more importantly, 
all the peaks fall into the same bin. This is explained con- 
sidering that, within this mass range, (i) stars are found 
to make the transition to the C-class at quite similar lu- 
minosities (i.e. Mboi ~ —4; see Fig. ||), and (ii) they span 
as carbon stars a luminosity interval which is compara- 
ble to the typical excursion of the low-luminosity dip, i.e. 
AMboi ~ 1 mag. 

In view of this similarity, we expect that the faint wing 
of the integrated CSLF is not significantly affected by pos- 
sible variations in the SFR occurring at ages older than 
~ 7 x 10 8 yr and, in particular, the location of the peak 
would remain unchanged. On the other hand, differences 
are evident among the distributions relative to more mas- 
sive (Mi > 2.5Mq) carbon stars. Hence, it follows that 
the high luminosity tail could reflect the details of the re- 
cent history of star formation, for ages comprised roughly 
within 7 x 10 8 and 10 8 yr. 



Fig. 11. Theoretical CSLFs for simple stellar populations. 
Each histogram corresponds to the expected distribution 
of carbon stars of the same age, i.e. evolved from progeni- 
tors with the same initial mass, as indicated. The shaded 
areas mark the position of the luminosity peak for each 
mass. Notice that all the histograms for M\ < 2.5M Q 
have the peak at the same luminosity bin. Calculations 
refer to TP-AGB models of initial metallicity Z = 0.008, 
for values of the dredge-up parameters, logTj drcd = 6.4 
and A = 0.5. Similar behaviours hold for other choices of 
the parameters. 



On the basis of the above considerations, let us now in- 
vestigate whether it is possible to reproduce the observed 
CSLF by only varying the adopted SFR, while keeping 
fixed all the other prescriptions. We start from the the- 
oretical distributions which are found not to fit the ob- 
served constraints with our standard prescription for the 
SFR [ip(t) = ipi, top panel of Fig. |l(J, for two different 
combinations of the dredge-up par ame ters (top and mid- 
dle panels of Fig. ||; see also Sects. 3d and |6.2|) . 

Concerning the (logT, drcd = 6.5, A = 0.65) case, shown 
in the top panel of Fig. ra, we can already conclude that 



the discrepancy cannot be removed, i.e. there is no way to 
shift the theoretical peak towards lower luminosities. In 
fact, as previously discussed, the faintest possible location 
of the peak is characteristic of all simple distributions of 
carbon stars with initial masses < 2.5M©, corresponding 
to ages older than ~7x 10 8 yr. This circumstance makes 
the lower limit of the peak luminosity almost invariant to 
any changes in the SFR over this age interval. We can 
then conclude that the luminosity location of the peak is 
indeed a strong calibrator of the efficiency parameter A. 
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Moreover, considering that the extension of the faint 
tail of the distribution, from the peak down to the end, 
is essentially controlled by the occurrence of the low- 



luminosity dip (see Sect. 5.4), it follows that it not possi- 
ble to populate the observed faintest bins, except by in- 
voking a larger extension of the subluminous stage (e.g. 
~ 1.75 mag) than assumed here (1.25 mag). 

Let us now analyse the (logT^ rcd = 6.4, A = 0.65) 
case, illustrated in the middle panel of Fig. |[ We aim at 
testing whether it is possible to calibrate a suitable law 
for the SFR, which is able to make the theoretical peak 
coincide with the observed one and to attain a general 
agreement of the overall features. After several trials, the 
best fit is obtained by approximating the SFR as the sum 
of a decreasing exponential and a gaussian function: 



tj}{t) oc exp ( — — | +- cxp 



t - t 

<7t 



(18) 



which is illustrated in the bottom panel of Fig. [ip(t) = 
ips]. The suitable value for the time-scale parameter is 
t s = 2.5 Gyr, producing quite a slow decline of the SFR to- 
wards younger ages. The effect of the slight decrease is not 
negligible, as it secures the reproduction of the faint end 
of the distribution, not reached otherwise with a strictly 
constant SFR. 

The gaussian function, centred at the age t = 5.5 x 
10 s yr and with a standard deviation cr t = 0.005, mimics 
the occurrence of a strong and short-lived burst around 
few 10 8 yr ago. This favours the contribution from car- 
bon stars belonging to relatively young stellar generations 
(Mi ~ 2.8Mq), so that the peak of the theoretical distri- 
bution is moved into the right luminosity bin. However, 
it is worth noticing that, although the overall luminosity 
range fits the observed one, an excess of luminous carbon 
stars (Mboj < —5.25 6) is predicted. To get a quan- 
titative estimate of corresponding error we also plot in 
Fig. [l2] the same theoretical distribution, but normalized 
to the total number of observed carbon stars brighter than 
Mbol = —5.25. In this case, the absolute agreement at high 
luminosities reflects into a considerable deficit of carbon 
stars fainter than Mboi > —5.25, covering a fraction of 
~ 37%, which is considerably smaller than the observed 
fraction, ~ 87%. The mismatch is really significant. 

Finally, let us consider the (logT b ircd = 6.4, A = 0.5) 
case, shown in the bottom panel of Fig. |8[ which is the 
best fit to the observed CSLF in the LMC. The theoretical 
distribution is obtained adopting a step- wise SFR [tp(t) = 
ipi in Fig. that corresponds to the occurrence of a 
dominant episode of star formation since the age % = 
5 x 10 9 yr, up to some more recent epoch t e . 

In practice, within the significant age range for the 
formation of carbon stars, ip(t) is constant up to f e , when 
the SFR is assumed to suddenly drop to very low (but 
not zero) values. From our analysis, it turns out that the 
best agreement with the observed high luminosity tail of 
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Fig. 12. Top panel: the theoretical LF of carbon stars 
(solid line) obtained with the dredge-up parameters as in- 
dicated, and a prescription for the SFR [ip(t) = 1P3] char- 
acterised by a recent burst at around 5.5 x 10 8 yr ago. 
The dotted line histogram represents the observed LF of 
carbon stars in the LMC. Bottom panel: The same as in 
the top panel, but with the theoretical histogram (dashed 
line) normalised to the number of observed carbon stars 
more luminous than Mboi = —5.25. 



the distribution (Afboi < —5.6) is obtained with t e ~ 5 x 
10 8 yr. Since the age defining the end of the burst is higher 
than that of the most massive (i.e. 4 — 5M©) carbon stars, 
their contribution remains quite modest, but sufficient to 
account for the bright end of the distribution at -Mboi ~ 
-6.8. 

If younger ages are assigned to i e , a slight excess of 
luminous carbon stars would show up. In this case, the 
discrepancy may be removed invoking, instead of a re- 
cent sharp decrease in the SFR, a more efficient envelope 
burning, or stronger mass-loss rates in massive TP-AGB 
stars. In this respect, it is worth noticing the high lumi- 
nosity bins of the distribution are crucially affected by 
the delicate interplay between convective dredge-up, en- 
velope burning, and mass-loss occurring in massive TP- 
AGB stars. Unfortunately, the present understanding of 
all these topics is incomplete so that theoretical predic- 
tions suffer from a degree of uncertainty. 

A further test is performed adopting a SFR law that 
is slowly decreasing - rather than constant - from an age 
of 5 Gyr up to now [if>(i] = ip2, middle panel of Fig. HG]: 



ip(t) oc exp(— t/t s ) 



(19) 
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Fig. 13. Top panel: the best fit (solid line) to the observed 
CSLF (dotted line) in the LMC. The theoretical distribu- 
tion is obtained with the adoption of the dredge-up pa- 
rameters as indicated, and following the prescription for 
the SFR (ip(t) = V>i) suggested by Bertelli et al. (1992). 
Bottom panel: The same as in the top panel, but assum- 
ing a SFR slowly decreasing with time [ip(t) = ip2] see also 
Fig. El. 

with a time-scale t s — 2.5 Gyr. The results, shown in 
Fig. [l^, are satisfactory for the faint wing of the distribu- 
tion, but both the peak and the high luminosity tail are 
somewhat less well reproduced than with the ip(t) — ipi 
prescription. 

From the above investigation it follows that the ob- 
served faint tail and the peak location of the CSLF cannot 
give stringent limits on the history of the SFR, as they 
are essentially determined by the dredge-up parameters. 
On the other hand, the bright wing of the distribution is 
sensitive to the details of the recent SFR. 

Let us now explore the sensitivity of the result to the 
IMF. In Fig. [l4] we present the resulting distributions of 
carbon stars adopting our best fit for the SFR (i.e. ipi 
shown in the top-panel of Fig. |To| ), while varying the law 
for the IMF according to three choices, namely: Salpeter 
(1955); Miller & Scalo (1979); and Ferrini et al. (1990). 
The first two prescriptions for the IMF are empirical, while 
the third one has a theoretical derivation. 

The differences among the predicted luminosity func- 
tions are indeed negligible. This can be understood con- 
sidering that: 

— the similarity of the predicted distributions for sim- 
ple stellar populations with turn-off masses 1 



Fig. 14. Predicted CSLFs according to different prescrip- 
tions for the IMF, as quoted in each panel (solid lines). 
The dotted line shows the observed distribution in the 
LMC. The adopted dredge- up parameters are (logT^ rcd = 
6.4, A = 0.5), and the SFR is fa. 

(M-JMq) <, 2.5 makes the results quite insensitive to 
the adopted IMF in this mass range, which indeed pro- 
vides the major contribution of carbon stars; 
— the small weight of the IMF at higher masses tends 
to cancel the differences in the results obtained with 
different prescriptions. 

In conclusion, it turns out that the observed distribu- 
tion of the carbon stars in the LMC is mainly determined 
by the the properties of the third dredge-up (mass-loss, and 
envelope burning), with the SFR and IMF producing a 
much weaker effect. 

7. The calibration of dredge-up parameters for 
the SMC 

The same kind of procedure described in the previous 
section is also applied to reproduce the CSLF in the 
SMC. In this case, we made use of synthetic calcula- 
tions for TP-AGB stars with initial chemical composition 
[Y = 0.240; Z = 0.004]. This values are consistent with 
present data for the abundances of young populations in 
the SMC. However, it is also found that the SMC clusters 
present a clear age-metallicity relation (Da Costa & Hat- 
dizimitriou 1998), with older clusters being systematically 
more metal poor. If we adopt the mean age-metallicity 
relation defined by the latter authors (i.e. excluding the 
apparently abnormal clusters L 113 and NGC 339), we 
find that it is necessary to go to ages as old as 5 Gyr 
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Fig. 15. Observed CSLFs (dotted line) in the LMC (top 
panel) and SMC (bottom panel), and the theoretical best 
fits (solid line). 



Concerning the broader extension in luminosity of the 
CSLF in the SMC if compared to the LMC, the following 
effects are remarked: 

— The faint end of the distribution is reached at lower 
luminosities in the case of lower metallicities mainly 
because of both an earlier onset of the third dredge-up, 
and the smaller number of dredge- up episodes required 
for the transition to the C-class (due to the lower orig- 
inal abundance of oxygen in the envelope). 

— As far as the bright end of the distribution is con- 
cerned, it is worth noticing that a good fit to the ob- 
served data for the SMC is obtained assuming a con- 
stant SFR over the entire significant age interval, with- 
out the recent drop (at an age of ~ 5 x 10 8 yr) invoked 
for the LMC. 

The CSLFs for both Magellanic Clouds are presented 
in the two panels of Fig. 



15. ft is remarkable that the 



present formulation for the third dredge-up allows a good 
fitting of the CSLFs of both Clouds without any dramatic 
change in the input parameters. In this context, the large 
difference between the CSLFs in the LMC and SMC are 
interpreted mainly as the result of a different metallicity. 
The different histories of SFR in both galaxies probably 
play only a minor role in their CSLFs. 



(corresponding to C stars with M\ ~ 1.25 M@) to find 
abundances reduced by a factor of 2. Therefore, it is pos- 
sible that Z = 0.002 evolutionary tracks would be more 
suitable to represent the oldest (and faintest) SMC car- 
bon stars. We remark, however, that the present investi- 
gation is intended to give a first hint about the behaviour 
of the dredge-up parameters at metallicities lower than 
Z = 0.008. Relevant transition luminosities as a function 
of the initial mass are displayed in Fig. |?], for two choices 
of the dredge- up parameters. 



The best fit is obtained adopting (log T[ 



drcd 



6.4, A 



0.65) in the simple case of a constant SFR and Salpeter's 
IMF as shown in the bottom panel of Fig.|l5|. It is inter- 
esting to notice that the proper value of is the same 
as for the LMC. This is in agreement with the indications 
from co mpl ete evolutionary calculations already discussed 
in Sect. 3^, that 7^ rcd is constant regardless of core mass 



and metallicity. 

Another aspect worthy of notice is that a greater ef- 
ficiency of the third dredge-up is required, if compared 
to the results for the LMC. This also agrees qualitatively 
with the indications from complete modelling of thermal 
pulses (e.g. Boothroyd & Sackmann 1988d) that a deeper 
envelope penetration is favoured at lower metallicities. As 
already mentioned, the effect of increasing A is that of 
shifting the peak location towards fainter luminosities. In- 
deed, this is one of the most evident differences between 
the distributions of carbon stars in the two Clouds. 



8. Conclusions 

In this paper we describe an improved treatment of 
dredge-up in synthetic TP-AGB models, based on the 
adoption of the parameter T^ lcd , i.e. the minimum base 
temperature required for dredge-up to occur. Every time 
a thermal pulse is expected during calculations, envelope 
integrations are performed to check whether the condition 
on the base temperature is satisfied. It follows that, con- 
trary to the test based on the constant M™ m parameter, 
with this new scheme the response depends not only on 
the core mass but also on the current physical conditions 
of the envelope (i.e. the surface luminosity peak, the ef- 
fective temperature, the mass, the chemical composition). 
Moreover, this method allows one to determine not only 
the onset, but also the possible shut-down of dredge-up 
occurring when the envelope mass is significantly reduced 
by stellar winds, without invoking further a priori assump- 
tions. 

This aspect is crucial. In fact, after the shut-down 
of the third dredge-up, surface and wind abundances of 
low-mass AGB stars are expected to freeze out till the 
end of the evolution, with important consequences for 
the chemical composition of planetary nebulae, which are 
ejected during the last evolutionary stages. Another point 
concerns the competition between envelope burning and 
dredge-up in the most massive TP-AGB stars. In partic- 
ular, if the shut-down of dredge-up occurs after envelope 
burning extinguishes, a star might undergo a late conver- 
sion to carbon star. The recent detection of a luminous 
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(Afbol ~ —6.8) dust-enshrouded carbon star in the LMC 
(van Loon ct al. 1998) seems to indicate such eventuality. 
This topic is the subject of a future study. 

The CSLFs in the Magellanic Clouds are investigated 
as possible calibrators of the dredge-up parameters, and 
possible indicators of the SFR history in the host galaxies. 
We briefly summarise here the basic points emerging from 
our analysis: 

— In principle, the observed CSLFs can provide indica- 
tions on the SFR history within the age interval from 
about 4 x 10 9 yr to ~ 10 8 yr. Nothing can be inferred 
for younger ages. 

— The faint end of the luminosity distribution is essen- 
tially determined by the temperature parameter X| dred . 

— The peak location depends only weakly on the SFR, 
and hence it is a stringent calibrator of the efficiency 
parameter A. 

— Although the peak location is essentially invariant with 
respect to the SFR, the shape of the carbon star dis- 
tribution may give hints on the history of the SFR in 
recent epochs (for ages between 7 x 10 8 and 10 8 yr 
ago). However, the bright wing is expected to be sen- 
sitive to the details of physical processes occurring in 
the most massive stars, i.e. envelope burning and mass 
loss. 

— The best fit to the observed CSLF in the LMC is ob- 
tained with A = 0.50, Tg led = 6.4, and a constant SFR 
up to an age of about 5 x 10 8 yr. This recent drop of 
the SFR is invoked to remove a slight excess of bright 
carbon stars otherwise predicted. 

— The best fit to the observed CSLF in the SMC is de- 
rived with A = 0.65, T b drod = 6.4, and a constant SFR 
over the entire significant age interval. 

— If compared with the calibration for the LMC, a 
greater efficiency of dredge-up is required to reproduce 
the fainter peak of the CSLF in the SMC. 
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